Weighing neutrinos in the presence of a running primordial spectral index 
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The three- year WMAP(WMAP3), combined with other cosmological observations from galaxy 
clustering and Type la Supernova (SNIa), prefers a non- vanishing running of the primordial spectral 
index independent of the low CMB multipoles. Motivated by this feature we study cosmological 
constraint on the neutrino mass, which severely depends on what prior we adopt for the spectral 
shape of primordial fluctuations, taking possible running into account. As a result we find a more 
stringent constraint on the sum of the three neutrino masses, m v < 0.76eV (2 a), compared with 
m v < 0.90eV (2 a) for the case where power-law prior is adopted to the primordial spectral shape. 



PACS number(s): 98.80.Cq 

The three year Wilkinson Microwave Anisotropy Probe 
observations (WMAP3) 0, S 0, fl 0] have marked an- 
other milestone on the precision cosmology of the Cosmic 
Microwave Background (CMB) Radiation. The simplest 
six-parameter power-law ACDM cosmology is in remark- 
able agreement with WMAP3 together with the large 
scale structure(LSS) of galaxy clustering as measured by 
2dF § and SDSS @ and with the Type la Supernova 
(SNIa) as measured by the Riess "gold" sample [8j and 
the first year SNLS [9|. This agreement between the 
above "canonical" cosmological model and observations 
can be used to test a number of possible new physics, 
such as the equation of state of dark energy, neutrino 
masses, time variation of fundamental constants, etc. 

Among them, the constraint on the neutrino or 
hot dark matter mass can be obtained from the free- 
streaming modification of the transfer function of the 
matter power spectrum. We should note, however, that if 
one allowed any shape of primordial spectrum, the free- 
streaming effect could easily be compensated by some 
nontrivial shape of the primordial spectrum, so that one 
cannot obtain sensible limit on the neutrino mass. That 
is, we can obtain a nontrivial bound on neutrino mass if 
and only if we adopt some prior on the shape of primor- 
dial power spectrum such as a simple power-law. From 
the above argument, we expect that as we allow more 
degrees of freedom on the primordial spectrum beyond a 
power-law, the constraint on the neutrino mass would be 
less stringent in general. 

As for the shape of the primordial power spectrum, 
it is noteworthy that a significant deviation has been ob- 
served by WMAP3 from the simplest Harrison-Zel'dovich 
spectrum, and that this feature is more eminent with the 
combination of all the currently available CMB, LSS and 
SNIa (dubbed the case of "All" in [5]). Moreover, a non- 
trivial negative running of the scalar spectral index a s , 
whose existence was studied even before WMAP epoch 
1C. Illl l , w as favored by the first-year WMAP papers 
12 . Il3l Il4j . But its preference was somehow diminished 
as corrections to the likelihood functions were made llal . 



However, the new WMAP3 data prefers again a negative 
running in the "All" combination |5|. 

If confirmed, a nonvanishing running of a s would 
not only constrain inflationary cosmology significantly 
[3, EI El IM iH, but also affect the cosmologi- 
cal constraint on the neutrino mass. In the LSS power 
spectrum, the effect of massive neutrino may be com- 
pensated by a nonvanishing running of primordial spec- 
trum. On the other hand, if negativeness of the run- 
ning is established, it will lead to even more stringent 
constraints on neutrino mass compared with fittings in 
the constant scalar spectral index(n s ) cosmology, because 
both of them lead to a damped power on small scales. 

The actual problem, however, cannot be solved by the 
above simple one-to-one correspondence, because the ef- 
fects of a nonvanishing neutrino mass on CMB is much 
less dramatic than a nonvanishing running a s . Alterna- 
tively, in the scales probed by CMB, especially near the 
third peak, there is a large de gen erac y b etween a s and 
the matter density Vt m 0, [HIM El, III]. 

Hence in the concordance analysis the correlation be- 
tween running and neutrino mass needs to be addressed 
in a combined study of CMB, LSS and SNIa, where 
SNIa helps significantly to determine the matter den- 
sity. We report the results of such a combined anal- 
ysis in this paper using Markov Chain Monte Carlo 
method in constraining the total neutrino mass, m„ = 
Y^i=i m vi = 94.4f2„/i 2 eV. Here f2„ is the density param- 
eter of the neutrino and h is the Hubble constant in unit 
of lOOkm/s/Mpc. 

To break possible degeneracy among the cosmological 
parameters, we make a global fit to the forementioned 
current data with the publicly available Markov Chain 
Monte Carlo package cosmomc [10Ll27|. Our most general 
parameter space is: 

p = (u> b , uj c , © s , t, m v , n s ,a s , log[10 10 A s ]) (1) 

where uib = Qbh 2 and w c = fl c h 2 are the physical baryon 
and cold dark matter densities relative to critical density, 
Og characterizes the ratio of the sound horizon and an- 
gular diameter distance, r is the optical depth and A s is 
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defined as the amplitude of initial power spectrum. The 
pivot scale for n s and a s is chosen at k = 0.05 Mpc -1 . 

Assuming a flat Universe and in terms of the Bayesian 
analysis, we vary the above 8 parameters and fit the the- 
ory to the observational data with the MCMC method. 
For CMB we have only adopted the WMAP3. The bias 
factors of LSS have been used as nuisance parameters 
and hence essentially we have used only the shapes of 
2dF and SDSS power spectra. As for the SNIa data, 
while the WMAP team uses both SNLS and Riess sam- 
ple simultaneously in their "All" dataset, here we adopt 
only one of them, namely, the Riess " gold" sample rather 
than combining with SNLS. This is because these two 
groups use somewhat different methods in their analysis 
and it would not be appropriate to put them together 
simply 28]. 

Regarding the first-year WMAP data Bridle et al. ^ 
found that the claimed preference of a negative a s was 
merely due to the lowest WMAP multipoles. In order 
to probe the sensitivity of the running to the lower mul- 
tipoles we analyze the running and neutrino properties 
using the CMB data with and without the contributions 
of lower multipoles which suffer from large cosmic vari- 
ance. Specifically, we truncate naturally at I = 24 given 
the current likelihood of WMAP 3 [3, [3, [3, . 

As a result we find a more stringent constraint on the 
neutrino mass in the presence of running compared with 
the analysis with constant n s . We also find that currently 
the preference of a negative running is fairly indepen- 
dent of the WMAP3 low CMB quadrupoles and hence 
relatively robust. 

In Table 1 we show the mean la (2a) constrains on the 
relevant cosmological parameters combining 2dF, SDSS 
and Riess "gold" sample with WMAP 3. We have ad- 
dressed the cases with/ without I < 24 CMB contribu- 
tions, with/without introducing a s and with/without 
massive neutrinos. The last column shows the cor- 
responding reduction of y 2 values compared with the 
power law ACDM cosmology. For normal LSS + SNIa 
+ WMAP in the 7 parameter fittings we get a s — 
— 0.0512lQ;Q4gQ at 2a. Our results are considerably less 
stringent than the "All" combination by WMAP team 
[B|, as we have not included SNLS and other CMB ob- 
servations. On the other hand the preference of non- 
vanishing a s is larger than 2a in both cases. The 
preference to negative running still remains even if we 
drop the WMAP3 I < 24 contributions, when we get 
a s = — 0.0717lQ Qggg at 2a. This may imply a negative 
running is indeed preferred nontrivially in the combina- 
tion of WMAP+LSS+SNIa, even without the presence 
of WMAP3 small / contributions. This can also be un- 
derstood through the Akaike information criterion (AIC 
[13]), which is defined as defined as 



AIC = -21n£ + 2k, 



(2) 



parameters of the model and N is the number of dat- 
apoints. From Table I we can find the reduction of 
A\ 2 = — 2Aln£ is larger than 2 in the case of running 
compared with the simple power law ACDM cosmology, 
hence we have smaller values of AIC with running, al- 
though currently not strong to decisive levels. 

The preference of a negative running can also be 
obtained in the two-dimensional posterior contours of 
n s — a s , as depicted in Fig. [TJ For the 7-parameter case 
with one additional parameter of a s , although the con- 
tour without WMAP3 low I contributions is larger than 
that in the left panel, a constant n s lies close to the 2a 
lines in both cases. The enlarged contour without I < 24 
is easily understood due to the n s — r degeneracy. 

The correlation between m„ and the shape of the pri- 
mordial spectrum is obvious, as can be seen from Table 
1 and Fig. [TJ In the presence of massive neutrinos the 
error bars on n s and a s get increased disregard with or 
without small I CMB contributions. And it is notewor- 
thy in the presence of massive neutrinos, a scale-invariant 
primordial spectrum is consistent with the observations 
at ~ 2(7 if we drop the small I WMAP3 contributions. 
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where C is the maximum likelihood, k is the number of 



FIG. 1: Two dimensional posterior constraints on the pri- 
mordial spectral index versus its running at k = 0.05 Mpc -1 , 
using 2dF, SDSS, SNIa and WMAP3 (with/ without I < 24 
CMB contributions) and with/without the presence of mas- 
sive neutrinos. 



We find that, while almost all of the remaining pa- 
rameters get less stringently constrained, the neutrino 
mass is an exception: a more tightened bound on rriy is 
achieved in the presence of a nonzero a s than the case 
with constant n s , as shown in Table 1 and in the one 
dimensional constraints as in Fig. O This has shown the 
fact that running is indeed strongly correlated with neu- 
trino mass, which is mainly due to the physics on LSS. 
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TABLE 1. Mean with la (2a) constrains on the spectral index, the running, and the neutrino mass based on LSS and SNIa 
with WMAP3(with/without I < 24 CMB contributions) and with/without introducing a running of the primordial spectral 
index a s and with/without massive neutrinos. The last column shows the corresponding reduction of % 2 values compared 

with the power law ACDM cosmology. 
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FIG. 2: One dimensional posterior constraints on the 
sum of neutrino masses, using 2dF, SDSS, SNIa and 
WMAP3(with/without / < 24 CMB contributions) and 
with/without introducing a running of the primordial spectral 
index a 3 . Note due to the non-gaussian distribution of m„, 
the LSS + SNIa + WMAP3 (black line) seems to indicate a 
nonzero neutrino mass, but it is not the case, as can be seen 
in Table 1. 
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FIG. 3: Two dimensional posterior constraints on the sum of 
neutrino masses versus matter density, using 2dF, SDSS, SNIa 
and WMAP3(with/without / < 24 CMB contributions) and 
with/without introducing a running of the primordial spectral 
index a s . 



In Fig. [3] we display the two dimensional posterior con- 
straints on the sum of neutrino masses versus matter 
density in the same case as Fig. [5] While the allowed 
parameter space for J7 m is significantly enlarged in the 
presence of running, neutrino mass is constrained strin- 
gently in cases with nonzero a s . The correlation between 
Q m and m v is rather strong and the accumulation of the 
observational data, such as SNAP, PLANCK and SDSS 
will help significantly to break the degeneracy, to detect 
the features of the primordial spectrum as well as the 
nature of neutrinos. 

It has been claimed that a running of the spectral in- 
dex will be excluded in the presence of SDSS Lyman a 
observations [3l|, [H, [33 , 34 1, the systematics of Lyman a 
data are relatively less constrained [l|, [23| and we leave 



this in a separate investigation to be reported in 35[ , and 
detailed analysis with other additional possible degenera- 
cies in [36j . 

For many years neutrinos have played a fundamen- 
tal role in both physics and astrophysics, and provided 
venues of new physics such as the parity violation and os- 
cillations with tiny masses. Surely neutrinos will continue 
to play a crucial role for our understanding of the Uni- 
verse. While the resulting reduction on neutrino mass 
in this paper is less than 0.2 eV, such an effect would 
hopefully help to change our understandings on the ulti- 
mate detection of neutrino mass with future cosmological 
surveys and the difference is already larger than the low 
limit of neutrino mass from oscillation experiments. On 
the other hand there are also some mild tensions in the 
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determinations of the background cosmological param- 
eters with current CMB, LSS and SNIa. We may still 
need some better understandings on each data-set before 
entering the precision cosmology 37J and in cases when 
all of the observations have similar tendencies in the pref- 
erence of a negative running one can hopefully get more 
eminent effects in the probe of neutrino mass. 

The distinctive feature probed in the current paper will 
also open new windows on relevant studies, such as prob- 
ing neutrino mass with a nonzero running in gravitational 
lensing surveys, N-body simulations in the presence of 
running and massive neutrinos. 
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